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The production of secondary electrons generated by carbon nanoparticles and pure water medium 
irradiated by fast protons is studied by means of model approaches and Monte Carlo simulations. It 
is demonstrated that due to a prominent collective response to an external field, the nanoparticles 
embedded in the medium enhance the yield of low-energy electrons. The maximal enhancement is 
observed for electrons in the energy range where plasmons, which are excited in the nanoparticles, 
play the dominant role. Electron yield from a solid carbon nanoparticle composed of fullerite, a 
crystalline form of C60 fullerene, is demonstrated to be several times higher than that from liquid 
water. Decay of plasmon excitations in carbon-based nanosystems thus represents a mechanism of 
increase of the low-energy electron yield, similar to the case of sensitizing metal nanoparticles. This 
observation gives a hint for investigation of novel types of sensitizers to be composed of metallic and 
organic parts. 


I. INTRODUCTION 

Radiotherapy is currently one of the most frequently 
used technologies to treat tumors, which are a major 
health concern [1]. However, this technique has a lim¬ 
itation which comes from the sensitivity of healthy tis¬ 
sues, surrounding the tumor, to radiation. To make the 
treatment more efficient, one needs to minimize the dose 
delivered to the healthy tissue, thus preventing harmful 
effects of radiation exposure. Therefore, approaches that 
enhance radiosensitivity within tumors relative to nor¬ 
mal tissues have the potential to become advantageous 
radiotherapies. A search for such approaches is within 
the scope of several ongoing multidisciplinary projects 
[2, 3], 

One of the most promising modern treatment tech¬ 
niques is ion-beam cancer therapy (IBCT) [4-6]. In this 
technique, radiation damage is initiated by fast ions in¬ 
cident on tissue. Propagating through the medium, the 
projectiles deposit their kinetic energy due to the ioniza¬ 
tion and excitation processes. Biodamage due to ionizing 
radiation involves a number of phenomena, which hap¬ 
pen on various spatial, time, and energy scales. The key 
phenomena can be described within the so-called multi¬ 
scale approach to the physics of radiation damage with 
ions (see reference [7] and references therein). As a re¬ 
sult of the interaction of projectiles with the medium, 
secondary particles, such as electrons, free radicals, etc., 
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are produced. By now, it is generally accepted that the 
vast portion of biodamage done by incident heavy ions is 
related to these secondary particles [7- 10]. Particularly, 
the low-energy electrons, having the kinetic energy from 
a few eV to several tens of eV, have been shown to act 
as important agents of biodamage [11, 12]. 

Metallic nanoparticles, especially those composed of 
noble metals, were proposed recently to act as sensitiz¬ 
ers in cancer treatments with ionizing radiation [13-17]. 
Such nanoagents delivered to the tumor region can boost 
the production of secondary electrons near the target 
[18, 19]. The enhanced production of low-energy elec¬ 
trons will also lead to an increase in the number of free 
radicals [20] as well as other reactive species, like hydro¬ 
gen peroxide H 2 O 2 [21], which can propagate from the 
cytoplasm to the cell nucleus. Thus, these species can 
deliver damaging impacts onto the DNA from the radi¬ 
ation induced damages associated with the presence of 
nanoparticles in other cell compartments [22]. An en¬ 
hanced production of the secondary species will lead to 
an increase of the relative biological effectiveness of ion¬ 
izing radiation. This quantity is defined as the ratio of 
the dose delivered by photons to that delivered by differ¬ 
ent radiation modalities, leading to the same biological 
effects, such as the probability of an irradiated cell inac¬ 
tivation. 

The physical mechanisms of enhancement of the elec¬ 
tron yield from sensitizing nanoparticles are still a de¬ 
bated issue. In the recent studies [23, 24], it was discov¬ 
ered that a significant increase in the number of emitted 
electrons due to irradiation of noble metal nanoparticles 
by fast ions comes from the two distinct types of collective 
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electron excitations. It was demonstrated that the yield 
of the 1 — 10 eV electrons is strongly enhanced due to the 
decay of plasmons, i.e. collective excitations of delocal¬ 
ized valence electrons in metallic nanoparticles. For elec¬ 
tron energies of about 10 — 30 eV, the dominating contri¬ 
bution to the electron yield arises from the atomic giant 
resonance associated with the excitation of d-electrons in 
individual atoms in a nanoparticle [23]. 

Excitation of plasmons by time-dependent external 
electric fields is a characteristic feature of not only metal¬ 
lic but also, to some extent, of carbon nanoscale systems. 
For instance, it is generally accepted that plasmon excita¬ 
tions dominate the spectra of photo- and electron impact 
ionization of fullerenes [25-30] and polycyclic aromatic 
hydrocarbons (PAHs) [31, 32]. 

In this paper, we demonstrate that the decay of plas¬ 
mons excited in carbon nanoparticles also plays a promi¬ 
nent role in the production of low-energy electrons. Due 
to the collective response to a time-dependent external 
electric field, these systems enhance the production of 
secondary electrons in a biological medium, in the en¬ 
ergy range where the plasmons play the dominant role. 
This is done by the calculation of spectra of secondary 
electrons ejected from a carbon nanoparticle composed 
of fullerite, a crystalline form of Ceo fullerene, irradiated 
by fast protons. The contribution of plasmon excitations 
to the electron production is evaluated by means of the 
plasmon resonance approximation [33-36]. The results of 
these calculations are compared to the model calculations 
based on the dielectric formalism [37] and Monte Carlo 
simulations [38, 39], carried out for pure water medium 
and for the medium with an embedded carbon nanopar¬ 
ticle. Utilizing and comparing different theoretical and 
numerical methods, we provide a recipe for evaluation of 
the electron production in the kinetic energy range from 
a few eV to thousands of eV. A single method does not al¬ 
low one to properly quantify the secondary electron yield 
in a broad energy range; thus, a combination of different 
approaches is required. 


II. THEORY AND COMPUTATIONAL DETAILS 
A. Plasmon resonance approximation 

The contribution of collective electron excitations to 
the ionization spectra of carbon nanoparticles is evalu¬ 
ated by means of the plasmon resonance approximation 
(PRA) (see references [33-36] and references therein). 
This approach postulates that the dominating contribu¬ 
tion to the ionization cross section in the vicinity of the 
plasmon resonance comes from collective electron excita¬ 
tions, while single-particle effects give a small contribu¬ 
tion compared to the collective modes [40, 41]. In the 
past, this approach has provided a clear physical expla¬ 
nation of the resonant-like structures in photoionization 
spectra [30, 34] and differential inelastic scattering cross 
sections [27, 28, 42, 43] of metallic clusters and carbon 


fullerenes irradiated by photons and fast electrons. 

To start with, we evaluate the plasmon contribution 
to the ionization spectrum of an isolated Cgo molecule. 
Within the utilized model, the fullerene is represented as 
a spherical ’’jellium” shell of a finite width, A R = R 2 — 
R\ , so the electron density is homogeneously distributed 
over the shell with thickness A R [26, 44, 45]. The chosen 
value, A R = 1.5 A, corresponds to the size of the carbon 
atom [44]. 

The interaction of a hollow system with a non-uniform 
electric field, created in collisions with charged projec¬ 
tiles, leads to the time-dependent variation of the elec¬ 
tron density appearing on the inner and outer surfaces 
of the hull as well as in its interior [36]. This varia¬ 
tion leads to the formation of a surface plasmon, which 
has two normal modes, the symmetric and antisymmetric 
[26, 44-46], and of a volume plasmon [41], which occurs 
due to a local compression of the electron density inside 
the shell. The detailed explanation of formation of dif¬ 
ferent plasmon modes can be found in references [34, 36]. 

The utilized approach relies on several parameters, 
which include the oscillator strength of the plasmon ex¬ 
citation, position of the plasmon resonance peak, and its 
width. The choice of these parameters can be justified 
by comparing the model-based spectra with either ex¬ 
perimental data or the results of more advanced ab initio 
calculations. As a benchmark of the utilized approach, 
the photo- and electron impact ionization cross sections 
of carbon-based systems, namely fullerenes and PAHs, 
were calculated recently [28, 30, 32, 43]. The results ob¬ 
tained for Ceo [28, 30, 43] agreed well with experimental 
data on photoionization [29] and electron inelastic scat¬ 
tering [28, 43]. Being a clear physical model which de¬ 
scribes collective electron excitations, the PRA has been 
proven to be a useful tool for interpretation of experi¬ 
mental results and making new numerical estimates. 

Within the PRA, the double differential inelastic scat¬ 
tering cross section of a fast projectile in collision with a 
hull-like system can be defined as a sum of three terms 
[28, 36] (hereafter, we use the atomic system of units, 
m e = \e\ = h = 1): 

d 2 (T p i dV s ) dV“) dV”) 
de 2 dflp 2 dc 2 dQ P2 d£ 2 dH P2 de 2 dH P2 

which describe the partial contribution of the surface (the 
two modes, s and a) and the volume (v) plasmons. Here 
£2 is the kinetic energy of the scattered projectile, p 2 its 
momentum, and fl p2 its solid angle. The cross section 
d 2 <Tpi/d£ 2 dH P2 can be written in terms of the energy loss 
Ae = £1 — £ 2 , of the incident particle of energy £ 1 . Inte¬ 
gration of d 2 er p i/dA£dn p2 over the solid angle leads to 
the single differential cross section: 

<?max 

= [ dn ^ = Jl. f d2gpl 

dAs J P2 dAed£} p2 P1P2 J q dAedfl P2 

Qmin 
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where pi is the initial momentum of the projectile and 
q = px — p 2 is the transferred momentum. Explicit ex¬ 
pressions for the contributions of the surface and vol¬ 
ume plasmons, entering equation (1), obtained within the 
plane-wave Born approximation, are presented in refer¬ 
ence [36]. The Born approximation is applicable since the 
considered collision velocities (v\ = 2 — 20 a.u.) substan¬ 
tially exceed the characteristic velocities of delocalized 
electrons in the fullcrene {v e « 0.7 a.u.). 

The surface and volume plasmon terms appearing on 
the right-hand side of equation (1) are constructed as a 
sum over different multipole contributions corresponding 
to different values of the angular momentum l: 

d V»> v cup r}° 

d£2d51p 2 1 (oj 2 — Ul^ 2 ) 2 + LJ 2 T^ 2 

J’.M v- ^r;-> 

d£ 2 df2 P2 “ \ ( 0,2 _ LUp'j 2 + u 2 ’ 

where i = s, a denotes the two modes of the surface plas¬ 
mon. Their frequencies are given by [36, 44]: 

^ = { lT 2lTl^ 1 + m + l)el+1 ) V2 % (4) 

where and ’+’ stand for symmetric (s) and antisym¬ 
metric (a) modes, respectively, and £ = R 1 /R 2 is the ra¬ 
tio of the inner to the outer radii of the shell. The volume 
plasmon frequency lo p , associated with the ground-state 
electron density po, is given by 


^ = \j r$(i - e) ’ <5) 

where N is the number of delocalized electrons involved 
in the collective excitation. In the case of a fullcrene 
C„, the number N of delocalized electrons represents the 
four 2s 2 2p 2 valence electrons from each carbon atom. 
Thus, we assume that 240 delocalized electrons of Ceo 
contribute to the formation of plasmons. 

In reference [42] it was shown that the excitations with 
large angular momenta have a single-particle rather than 
a collective nature. With increasing /, the wavelength of 
plasmon excitation, A p i = 2i tR/1, becomes smaller than 
the characteristic wavelength of the delocalized electrons 
in the system, A e = 2i:/y/2e. Here e is the characteristic 
electron excitation energy in the cluster, e ~ J p , and I p is 
the ionization threshold of the system (I p (Cqo) ~ 7.5 eV 
[25]). In the case of the Cgo fullcrene, the estimates 
show that the excitations with l > 3 are formed by sin¬ 
gle electron transitions rather than by the collective ones. 
Therefore, only terms corresponding to the dipole {I = 1), 
quadrupole (/ = 2 ) and octupole (7 = 3) plasmon terms 
have been accounted for in the sum over l in equation (3). 

Following the methodology utilized in reference [43], 
we assume that the ratio 7 ; = T; /uji of the width of the 


TABLE I. Peak positions of the surface and the volume plas¬ 
mon modes as well as their widths used in the present calcu¬ 
lations. All values are given in eV. 



1 = 1 

l = 2 

l = 3 


19.0 

25.5 

30.5 

p( s ) 

1 1 

11.4 

15.3 

18.3 


33.2 

31.0 

29.5 

p(a) 

1 l 

33.2 

31.0 

29.5 

lOp 


37.1 


r»(v) 

1 l 


37.1 



plasmon resonance to its frequency equals to 7 ^ = 0.6 
for all multipole terms of the symmetric mode, and to 
= 1.0 for the antisymmetric mode. These values 
were utilized previously to describe experimental data 
on photoionization [29] and electron inelastic scatter¬ 
ing [28, 43] of gas-phase Ceo- The value 7 ^ = 0.6 
is also close to the numbers obtained from the earlier 
photoionization and electron energy loss experiments on 
neutral Ceo [25, 27]. The value 7 ^ = 1.0 is consis¬ 
tent with the widths of the second plasmon resonance 
observed in the photoionization of CgJ (q = 1 — 3) ions 
[47]. For the volume plasmon, we consider the ratio 

7 ; 1 '* = T^/wp = 1.0. The values of the plasmon res¬ 
onance peaks and the widths are summarized in Table I. 


B. Dielectric formalism 

The secondary electron production in a pure water 
medium as well as in a carbon nanoparticle was investi¬ 
gated by means of a model approach based on the dielec¬ 
tric formalism [37]. This method relies on experimental 
measurements of the energy-loss function of the target 
medium, Im[—l/e(w, q)], where e(cu, q) is the complex di¬ 
electric function, with uj and q being the energy and the 
momentum transferred to the electronic excitation, re¬ 
spectively. In reference [48], this approach was used to 
obtain spectra of secondary electrons generated in liq¬ 
uid water by energetic ions. An alternative method to 
calculate the impact ionization cross sections of various 
biological media was proposed recently [49, 50]. Instead 
of calculating the exact energy-loss function and ioniza¬ 
tion threshold for different electronic shells of a molecule 
composing the target medium, this approach aims at cal¬ 
culating the mean value of the binding energies for several 
outer shells. It is assumed that ionization of these shells 
happens if the energy transferred to the medium exceeds 
this mean value of the binding energies [49]. The formal¬ 
ism presented allows one to calculate the cross sections 
not only for liquid water but also for a real biological 
medium containing sugars, amino acids, etc. In particu¬ 
lar, it was utilized recently [51] to study ionization and 
energy deposition in different subcellular compartments, 
such as cell nucleus and cytoplasm, due to proton irra- 
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diation. In this work, we apply this formalism to study 
the electron production from a nanoparticle composed of 
fullerite. 


C. Monte Carlo simulation of secondary electron 
yield 

Monte Carlo simulations of secondary electron produc¬ 
tion in a nanoparticle were performed using Geant4, ver¬ 
sion 9.6 patch 1 [38, 39]. The simulation geometry con¬ 
sisted of a 50 nm diameter spherical nanoparticle of vari¬ 
able material placed at the center of a 5 /rm world of 
liquid water. A 4 /rm sided cube was included to allow 
the use of different secondary particle production thresh¬ 
olds in different regions in order to optimize execution 
times. Monoenergetic protons propagating from a point 
source were incident from the edge of the nanoparticle. 

The material of the nanoparticle was simulated as liq¬ 
uid water or a customized fullcrene material alternatively. 
The fullerene material properties were set by scaling the 
density of the Geant4 element carbon according to the 
calculated density of a face-centered cubic (fee) structure 
of fullerite. 

The Low Energy Electromagnetic Physics Package 
[52], using the Livermore Data Libraries, was selected 
to model the interactions of electrons and photons in 
the nanoparticle. Models describing proton interac¬ 
tions in the nanoparticle were selected following the 
Geant4 advanced example ’’Microdosimetry”. The ion¬ 
ization model implemented for protons was the Geant4 
” BragglonGas” model, valid for protons kinetic en¬ 
ergy up to 2 MeV, while the Bethe-Bloch model was 
adopted for higher energies. In the nanoparticle, nu¬ 
clear stopping power was modeled using the Geant4 
”ICRU49NucStopping” model. The multiple scattering 
was modeled for all charged particles with the Geant4 
”UrbanMsc95” model [52]. Atomic de-excitation (fluo¬ 
rescence and Auger electrons) was modeled as well [53]. 
Secondary electron production from the nanoparticle is 
limited to the electrons with kinetic energy greater than 
250 eV as this is the low-energy limit of validity of the 
Livermore Data Libraries [54]. 

The Geant4-DNA Very Low Energy extensions [55] 
were adopted in liquid water surrounding the nanoparti¬ 
cle to model in detail particle interactions down to a few 
eV scale. Physical interactions modeled for protons in the 
water sphere were G4DNAExcitation, G4DNAIonisation, 
and G4DNAChargeDecrease. The models used are the 
default Geant4-DNA model classes. 

The simulations in this study modeled the interactions 
of 1 MeV protons generated from one position and in one 
direction incident on the nanoparticle. Secondary elec¬ 
trons were produced in the nanoparticle with a cut of 250 
eV. The cut is the threshold of production of secondary 
particles. Below the cut, secondary electrons are not pro¬ 
duced and their energy is deposited locally, while above 
the cut, secondary electrons are produced and tracked in 


the nanoparticle and in the surrounding medium. The ki¬ 
netic energy spectra of secondary electrons escaping the 
nanoparticle were retrieved and the spectra at creation 
was compared directly to the same physical quantity cal¬ 
culated by means of the analytical model. The propor¬ 
tion of escaping secondary electrons produced within the 
fullerite-like nanoparticle was 98.5%. 

III. RESULTS AND DISCUSSION 

A. Electron production by an isolated C6o molecule 
due to the plasmon excitation mechanism 



FIG. 1. Upper panel: contribution dtjpi/dAg: of the plas¬ 
mon excitations to the single differential cross section of Ceo 
fullerene irradiated by fast protons of different incident ener¬ 
gies as a function of the energy loss. Lower panel illustrates 
the contribution of different plasmon excitations to the cross 
section dcrpi/dAe of C6o irradiated by a 1 MeV proton. 

The upper panel of Figure 1 illustrates the single dif¬ 
ferential cross section dcr p i/dAe calculated by means of 
the PRA for the Ceo fullerene irradiated by fast pro¬ 
tons of different incident energies as indicated. The pre¬ 
sented spectra comprise contributions of both the sur¬ 
face and volume plasmon excitations of different angular 
momenta l. As mentioned in Section II A, we have ac¬ 
counted for the dipole (l = 1), quadrupole (l = 2), and 
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FIG. 2. Relative contribution of different multipole terms to the single differential cross section d<7 p i/dAe of Ceo fullerene 
irradiated by 0.1, 1, 10 MeV protons as a function of the energy loss. 


octupole (Z = 3) plasmon terms because the excitations 
with higher angular momentum are formed by single elec¬ 
tron transitions rather than by the collective ones. The 
contribution of different plasmon modes to the spectrum 
of Cgo irradiated by a 1 MeV proton is illustrated in the 
lower panel of Figure 1. The main contribution to the 
cross section comes from the symmetric mode of the sur¬ 
face plasmon, whose relative contribution exceeds that 
of the volume plasmon by about an order of magnitude. 
The similar trend was observed recently studying elec¬ 
tron production by noble metal nanoparticles [23, 24]. 
Thus, the leading mechanism of electron production by 
sensitizing nanoparticles due to the plasmon excitations 
should be related to the surface term but not to the vol¬ 
ume one. 

Figure 1 demonstrates that the amplitude and the 
shape of the plasmon resonance depend strongly on the 
kinetic energy of protons. It was shown previously [42] 
that the relative contributions of the quadrupole and 
higher multipole terms to the cross section decrease sig¬ 
nificantly with an increase of the collision velocity. At 
high velocities, the dipole term dominates over the con¬ 
tributions of larger l, since the dipole potential decreases 
slower at large distances than the higher multipole po¬ 
tentials. To illustrate this effect, we have plotted the 
partial contributions of different multipole modes which 
are excited due to irradiation by 0.1, 1, and 10 MeV pro¬ 
tons. These dependencies are presented in Figure 2. For 
the sake of clarity, the cross sections, which represent the 
sum of three multipole contributions, have been normal¬ 
ized to unity at the point of maximum. Thus, one can 
compare directly the relative contribution of the different 
terms to the cross section at different incident energies. 
A prominent interplay of the different multipole terms at 
the lowest incident energy (left panel) results in a shift 
in the position of the maximum of the cross section. 

To quantify the production of electrons in collision 
with a nanoparticle, we redefine the cross section du/dAe 
as a function of the kinetic energy W of emitted elec¬ 
trons. This quantity is related to the energy loss via 
W = Ae — I p , where I p is the ionization threshold of the 
system. The first ionization potential of the Cgo fullerene 


approximately equals to 7.5 eV [25]. 



FIG. 3. Single differential cross section da/AW of the Ceo 
fullerene (thick solid and dash-dotted black curves) and of 
a water molecule (thin solid and dashed blue curves) irradi¬ 
ated by a 1 MeV proton as a function of the kinetic energy 
of emitted electrons. Thick solid (black) curve illustrates the 
contribution of the plasmon excitations to the emission spec¬ 
trum from Ceo • Thin solid and dashed (blue) curves represent 
the results obtained within the dielectric formalism by Scifoni 
et al. [48] and de Vera et al. [49], respectively. Symbols repre¬ 
sent the cross section of a single C atom calculated by means 
of BEA, multiplied by 60. 

Figure 3 shows the cross section dcr/dlF of Ceo (thick 
solid and dash-dotted black curves) and of a water 
molecule (thin solid and dashed blue curves) irradiated 
by a 1 MeV proton as a function of the kinetic en¬ 
ergy of emitted electrons. The results for water ob¬ 
tained within the dielectric formalism are taken from ref¬ 
erences [48, 49]. The thick solid curve demonstrates the 
contribution of the plasmon excitations to the spectrum 
of Ceo, dcrpi/d W, calculated within the PRA approach. 
The dash-dotted curve represents the results obtained 
within the dielectric formalism. In the latter case, we 
took the experimental optical energy-loss function of ful- 
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lerite [56] and calculated the mean binding energy of the 
outer-shell electrons. The binding energies of the valence 
orbitals of Cgo were taken from the ab initio calculations 
of Deutsch et al. [57]. Symbols show the cross section 
da/dW for the 1 MeV proton impact of a single car¬ 
bon atom calculated by means of the binary encounter 
approximation (BEA) [58, 59], multiplied by 60. The re¬ 
sults of the calculations based on the dielectric formalism 
agree well with those within the BEA at the energy of 
about 20 cV and above. This indicates that the emission 
of electrons with kinetic energy of about several tens of 
eV takes place via single-electron excitations of the sys¬ 
tem. The plasmon excitations dominate the spectrum at 
lower energies, i.e. in the vicinity of the plasmon reso¬ 
nance, while this contribution drops off at higher ener¬ 
gies of emitted electrons. In the energy range where the 
plasmons are excited, single-particle effects give a small 
contribution as compared to the collective modes. At 
higher energies, the collective excitation decays to the in¬ 
coherent sum of single-electron excitations. Note that at 
lower electron energies (from 1 to approximately 20 eV) 
the BEA-based results start to deviate significantly from 
that of the dielectric formalism. This deviation indicates 
that the BEA is not applicable for the description of 
low-energy electron emission, since these electrons are 
produced in distant rather than in binary collisions. In 
this energy range, the PRA approach better describes 
the low-energy electron emission since it accounts for the 
collective electron effects omitted in other models. 


B. Electron production by a large carbon 
nanoparticle 


In the previous section, we have calculated the sin¬ 
gle differential cross section for an isolated Cgo molecule 
within the PRA approach and the dielectric formalism. 
Now, we we apply these methods as well as the Monte 
Carlo scheme to study the production of secondary elec¬ 
trons by a large solid carbon nanoparticle whose density 
corresponds to that of fullerite, the crystalline form of 

C 60 - 

The single differential cross section da/dW can be re¬ 
lated to the probability to produce N secondary electrons 
with kinetic energy W, in the interval dW, emitted from 
a segment Ax of the trajectory of a single ion [7, 60]: 


d N(W) 
dW 


= nAx 


da 

dW 


( 6 ) 


where n is the atomic density of a system of compounds, 



200 400 600 800 1000 


W (eV) 

FIG. 4. Number of electrons per unit energy produced by 
irradiation of a 50 nm carbon nanoparticle by a single 1 MeV 
proton (black curves and filled circles). Blue curves represent 
the number of electron generated in the equivalent volume 
of liquid water. Solid and dashed blue curves represent the 
results obtained within the dielectric formalism by Scifoni et 
al. [48] and de Vera et al. [49], respectively. Open circles 
illustrate this quantity obtained on the basis of Monte-Carlo 
simulations using the Geant4-DNA simulation tool. 


(i) C6o molecules in fullerite are packed in the fee crys¬ 
talline lattice, and (ii) a unit cell is composed of four Ceo 
molecules. Knowing the lattice parameter of fullerite, 
a = 1.417 nm, and the mass of a single carbon atom, 
me = 12 u, we have calculated the density of the ful¬ 
lerite crystal: 

p(fullerite) = 60 m C _ ^ gg g/ern 3 . (8) 

a 6 

Utilizing these values, we have obtained the atomic den¬ 
sity of fullerite: 


P 

At at m at 


( 7 ) 


n(fullerite) = AAli£Ut£> 
60 • me 


1.4- 10 21 cm" 3 , 


(9) 


with p being the mass density of a target, _/V at the number 
of atoms in the target compound, and m a t the atomic 
mass. 

As a case study, we have considered a nanoparticle of 
50 nm in diameter. In the calculations, we assumed that 


which is by about an order of magnitude smaller than 
that of water, n(water) = 3.34 • 10 22 cm -3 . 

In Figure 4, we compare the electron yield from a 
50 nm spherical carbon nanoparticle and from the equiv¬ 
alent volume of pure water medium. We have calculated 
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the number of electrons per unit energy produced due to 
irradiation by a 1 MeV proton. Thick black curve rep¬ 
resents the contribution of collective electron excitations 
estimated by means of the PRA. The dash-dotted black 
curve shows the number of electrons estimated by means 
of the dielectric formalism. Filled and open symbols rep¬ 
resent the results of the Monte Carlo simulations carried 
out by means of the Geant4 tool for the carbon nanopar¬ 
ticle and pure water medium, respectively. Thin solid 
and dashed blue curves represent the results of recent 
calculation for liquid water obtained within the dielec¬ 
tric formalism [48, 49]. Note that in the Monte Carlo 
simulations, we did not simulate the crystalline lattice of 
fullerite explicitly but the material properties were set by 
scaling the density of the Geant4 element carbon accord¬ 
ing to the calculated density p (fullerite). 

Comparative analysis of the spectra at low kinetic en¬ 
ergy of emitted electrons (the upper panel of Figure 4) 
demonstrates that the number of electrons with the en¬ 
ergy of about 10 eV, produced by the carbon nanopar¬ 
ticle via the plasmon excitation mechanism , is several 
times higher than that created in pure water. The en¬ 
hancement of the yield of low-energy electrons may in¬ 
crease the probability of the tumor cell killing due to the 
double- or multiple strand break of the DNA [7]. Simi¬ 
lar to the case of noble metal nanoparticles [13-17], the 
use of carbon-based nanostructures in cancer treatments 
with ionizing radiation can thus produce the sensitization 
effect. As it was shown recently [23, 24], the number of 
electrons with the energy of about a few eV produced by 
the noble metal (gold and platinum) nanoparticles via 
the plasmon excitation mechanism exceeds that gener¬ 
ated in the same volume of liquid water by an order of 
magnitude. In the case of a carbon nanoparticle, the elec¬ 
tron yield reaches the maximum at higher electron ener¬ 
gies, namely at about 10 eV. Assuming this, one can con¬ 
sider novel metal-organic sensitizing nanoparticles, where 
collective excitations will arise in both parts of the sys¬ 
tem. A proper choice of the constituents will allow one 
to tune the position of the resonance peaks in the ioniza¬ 
tion spectra of such systems and, subsequently, to cover 
a broader kinetic energy spectrum of electrons emitted 
from such nanoparticles. The fabrication of new, more 
efficient types of sensitizers would allow one to signifi¬ 
cantly advance modern techniques of cancer treatment 
with ionizing radiation. 

In the case of electrons with higher kinetic energy (the 
lower panel of Figure 4), the effect done by the carbon 
nanoparticle (filled symbols and dash-dotted black curve) 
is also more prominent as compared to pure water (open 
symbols and dashed blue curve), as follows from both 
the calculations based on the dielectric formalism and 
the Monte Carlo simulations. As discussed above, the 
contribution of the plasmon excitations rapidly decreases 
at the energies exceeding approximately 30 eV. The PRA 
accounts only for collective electron excitations that dom¬ 
inate the ionization spectra at low energies. At higher 
energies, the plasmons decay into the incoherent sum of 


single-electron excitations whose contribution is the most 
prominent in this energy region. 



FIG. 5. Yield enhancement from the 50 nm carbon nanopar¬ 
ticle as compared to pure water medium. Solid and dashed 
blue lines show the enhancement due to the plasmon excita¬ 
tions as compared to the results obtained within the dielectric 
formalism by Scifoni et al. [48] and de Vera et al. [49], respec¬ 
tively. Open symbols illustrate the plasmon-based enhance¬ 
ment compared to the results of Monte Carlo simulations. 
The enhancement estimated solely by means of the dielectric 
formalism and the Monte Carlo simulations in a broader ki¬ 
netic energy range is shown in the inset by the dash-dotted 
curve and filled symbols, respectively. 

In order to quantify the difference in electron produc¬ 
tion by the carbon nanoparticle and by an equivalent 
volume of pure water, we have calculated the relative 
enhancement of the electron yield from the nanoparti¬ 
cle as compared to water. This quantity is presented 
in Figure 5. The main figure shows the enhancement 
which was calculated by comparing the contribution of 
the plasmon excitations, obtained within the PRA, to 
the electron yield from pure water calculated by means 
of the dielectric formalism (solid and dashed blue curves) 
and Monte Carlo simulations (open symbols). Depend¬ 
ing on the data to be chosen as a reference, the collective 
electron excitations result in 2 to 3 times greater number 
of emitted 10 eV electrons as compared to the case of 
water. This effect is less pronounced than the enhance¬ 
ment done by small noble metal nanoparticles which can 
produce up to 15-20 times greater number of electrons 
via the plasmon decay mechanism as compared to water 
[23, 24]. On the other hand, this enhancement results 
in an excessive emission of the very low-energy electrons 
of about a few eV, while the carbon-based nanoparticle 
can enhance the yield of more energetic electrons. For 
the sake of completeness, we also demonstrate the en¬ 
hancement done by the carbon nanoparticle in a broader 
kinetic energy range (see the inset of Figure 5). For that 
purpose, we have compared the electron yields from the 
two systems calculated by means of the dielectric formal¬ 
ism (dash-dotted curve) and also from the Monte Carlo 








simulation (filled symbols). The two approaches lead to 
a similar result, namely that the carbon nanoparticle en¬ 
hances the number of energetic (of about hundreds of eV 
up to 1 keV) secondary electrons by about 50%. 

The analysis performed demonstrates that a single the¬ 
oretical or numerical approach does not allow one to 
properly quantify the secondary electron yield in a broad 
kinetic energy range, from a few eV up to a few keV. 
Thus, one needs to utilize a combination of different 
methods to achieve this goal. The calculated spectra of 
secondary electrons can further be used as the input data 
for investigation of radiobiological effects by means of the 
multiscale approach to the physics of radiation damage 
with ions [7]. This approach has the goal of developing 
knowledge about biodamage at the nanoscale and molec¬ 
ular level and finding the relation between the charac¬ 
teristics of incident particles and the resultant biological 
damage. 

IV. CONCLUSION 

We have analyzed numerically the production of elec¬ 
trons by carbon nanoparticles irradiated by fast protons. 
The study has been carried out by means of the model 
approaches based on the plasmon resonance approxima¬ 
tion and the dielectric formalism, as well as by means 
of Monte Carlo simulations. It has been demonstrated 
that due to the prominent collective response to a time- 
dependent external electric field, carbon-based nanopar¬ 
ticles enhance the production of low-energy electrons via 


the plasmon excitation mechanism. 

The contribution of plasmons to the electron produc¬ 
tion from a carbon nanoparticle has been compared to 
the results of model calculations, based on the dielec¬ 
tric formalism, as well as to the results of Monte Carlo 
simulations for pure water medium. It has been shown 
that the number of the low-energy electrons (with the 
kinetic energy of about 10 eV) produced by a 50 nm car¬ 
bon nanoparticle is several times higher than that emit¬ 
ted from pure water. Similar to the case of sensitizing 
metallic nanoparticles, the decay of the plasmon excita¬ 
tions formed in carbon nanostructures represents an im¬ 
portant mechanism of generation of low-energy electrons. 
This observation gives an opportunity to fabricate new 
types of sensitizers, composed of the metallic and the or¬ 
ganic parts, where the plasmon excitations will arise in 
both parts of the system. As a result, it will become pos¬ 
sible to cover a broader kinetic energy range of electrons 
emitted from such systems, as compared to currently pro¬ 
posed nanoagents, and, subsequently, to improve modern 
techniques of cancer treatment with ionizing radiation. 
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